Role of Spin Diffusion in Current-Induced Domain Wall Motion 
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Current-induced spin torque and magnetization dynamics in the presence of spin diffusion in 
magnetic textures is studied theoretically. We uncover an additional torque on the form ~ V 2 [M x 
(u- V)M], where M is the local magnetization and u is the direction of injected current. This torque 
is inversely proportional to the square of the domain wall width (w and strongly depends on 
the domain wall structure. Whereas its influence remains moderate for transverse domain walls, it 
can significantly increase the transverse velocity of vortex cores. Consequently, the spin diffusion 
can dramatically enhance the non-adiabaticity of vortex walls. 
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The electrical control of the magnetic state of 
nanoscale heterostructures ll( such as magnetic domain 
walls [H and vortex cores Q is attracting increasing in- 
terest as a promising mechanism for innovative memory 
devices 0|. Identifying the nature of the torque exerted 
by the injected current on the domain wall itself has con- 
stituted a stimulating challenge resultingin the obser- 
vation of unique dynamical behaviors |2t6| and raising 
seminal questions concerning the transport of itinerant 
spins in inhomogeneous magnetic textures [l"12|. The 
most widely accepted form of the spin transfer torque ex- 
erted by a charge current on a magnetic texture M(r, t) 
is @,| 

T = 6j(u- V)M-^— Mx (u- V)M (1) 

M s 

where bj is the adiabatic spin torque, /3 describes the non- 
adiabaticity of the spin torque and u is the direction of 
current injection. The non-adiabaticity /? is generated by 
different mechanisms such as spin relaxation 0, EI| and 
magnetic texture-induced spin mistracking [9Hll| . In ad- 
dition to non-adiabaticity, it has recently been found that 
the Gilbert dampin g a can also be affected (enhanced) 
by the spin texture 12j. From the viewpoint of domain 



wall dynamics, the longitudinal (transverse) velocity of 
transverse domain walls (vortex cores) is controlled by 
the ratio of (3 to a @- Therefore, experimental efforts 
have been expended in accurately determining /3 and a 
for a wide range of magnetic materials and domain wall 
widths 0-[|. 

Recent experiments have shown that this ratio depends 
on the domain wall structure (Bloch, Neel or Vortex 
wall). Thomas et al. |f|, Eltschka et al. [l3[ and Hcyne 
et al. (l4j have found that vortex cores exhibit a much 
larger non-adiabaticity (/? w 8a to 10a) compared to 
transverse domain walls (/? w a). The authors attribute 
these large non-adiabaticities to the narrow character of 
domain wall width in vortex cores (~10 nm). As a mat- 
ter of fact, Tatara et al. and Xiao et al. Q demon- 



strated that in sharp domain walls the itinerant spin can- 
not adiabatically follow the local spin texture, resulting 
in an enhancement of the non-adiabaticity. On the other 
hand, Burrowes et al. fl5l ] have tested a very sharp trans- 
verse wall of about 1 nm using FcPt nanowircs and found 
that such a narrow domain wall docs not cause a signifi- 
cant increase in the non-adiabaticity, (3 « a. Theoretical 
investigations by Bohlens and Pfannkuche recently 
showed that the non-adiabatic torque has damped oscil- 
latory behavior when increasing the domain wall width, 
which may account for the non-adiabaticity enhancement 
for sharp domain wall widths. However, this model is ap- 
plied to Bloch walls only and does not readily explain the 
observed differences between transverse walls and vortex 
cores. 



The experimental observations SEMI indicate that 
the nature of the non-adiabatic spin torque exerted on 
sharp magnetization patterns is related to their dimen- 
sionality. Indeed, whereas a magnetic transverse wall 
varies along one direction only (c^M ^ 0, <9 y M = 0), 
the magnetization of a vortex core varies along two di- 
rections (dxM ^ 0, <9 y M ^ 0). Therefore, one ap- 
proach to explain the different non-adiabaticities of sharp 
transverse domain walls and vortex cores is to consider 
a mechanism that couples both a; and y directions. We 
have recently shown that a transverse spin current caused 
by anomalous Hall effect increases the transverse veloc- 
ity of isolated vortex cores while lea ving the transverse 
domain walls essentially unchanged [16J. However, its 
contribution is of the order of the damping constant a, 
which is insufficient to explain the results observed in 
Refs. @, [H Ell- 
in this letter, we demonstrate that spin diffusion gives 
rise to an additional spin torque that contributes to the 
current-driven velocity in the case of both Bloch walls 
and isolated vortex cores. However, whereas this addi- 
tional component affects the velocities of Bloch walls only 
moderately, it can dramatically increase the velocities of 
vortex cores. The itinerant electrons evolving in mag- 



2 



netic textures are described by the one-electron Hamil- 
tonian 



H 



2m 



— <x-M(r,i), 



(2) 



where p is the momentum, a is the Pauli matrix, M(r, t) 
is the time- and spatial-dependent magnetization (|M| = 
M s ) and J cx is the exchange coupling energy. Using 
Ehrenfest relation dtm = dt (cr) = i/h([&, H}), we ob- 
tain the spin continuity equation for the spin density m 

d t m=-V-J l — <5mxM-— , (3) 

where J = (v (g> <r) is the spin current tensor, r ox = 
h/2J ex is the spin precession time 10~ 15 s to 10~ 14 
s), r s f is the phenomenological spin-flip relaxation time 
(w KT 13 s to 1CT 12 s). We assume m = f^M + 5m, 
n s (5m) being the (non-) equilibrium spin density (n s ~ 
10" 2 M S ). 

Notice that Eq. ([3]) should in principle include spin 
dephasing contribution, arising from the destructive in- 
terference of non-equilibrium spins with different wave 
vector direction. Microscopic investigations using real- 
istic Fermi surfaces have shown that this effect destroys 
the transverse component of itinerant spins within a few 
monolayers in strong ferromagnets 17J. Disregarding 
this mechanism renders the algebra more tractable with- 
out qualitatively modifying the conclusion of the present 
Letter. 

As mentioned above, a promising candidate to explain 
the experimental observations is a mechanism that would 
link both x-y directions. In diffusive systems such as 
ferromagnetic metallic nanowires, the local momentum 
scattering tends to counteract the influence of the ex- 
ternal electric field resulting in an additional diffusion 
term in both charge and spin currents [181 ] that couples 
x and y directions. Moreover, in order to track possi- 
ble influence of sharp magnetic textures on the damping 
constant a, we also incorporate the correction proposed 
by Zhang and Zhang [12j. These authors have shown 
that in a time- and spatial-dependent magnetic texture, 
local spin pumping (or spin-motive force - SMF) results 
in a spatial-dependent tensor form of the damping [l2j ]. 
In the drift-diffusion approximation, in the presence of a 
local spin motive force (see Ref. 12| for details), the i-th 
component of the spin current reads 

Ji = -b.,uM - D diSm + Yrd t M x fcM. (4) 

Here bj = hbPGqE / 'eM s and 77 = gfi B ^G /4e 2 M s , P is 
the spin polarization, hb is the Bohr magnetron, Go the 
electrical conductivity, Dq is the diffusion coefficient and 
Ui is the i-th direction of the injected current. In Eq. (j4}, 
the first term is the adiabatic spin current driven by the 
current density j = GqEu, the second term arises from 
carrier diffusion and the last term arises from the SMF 



[12(1 . Notice that spin diffusion itself reduces the mag- 
nitude of SMF-damping dependin g on the spin diffusion 
length and the domain wall width 1!| 2(|. However, the 
analytical treatment of this effect is complex, so in the 
present work we assume a moderate value of the SMF 
parameter 77 to qualitatively account for this reduction. 

By injecting Eq. ((4]) into Eq. ([3]), we obtain the diffu- 
sion equation for the non-equilibrium spin density <5m 

1 r 8m 
5m x M H 

X M S T s f 

bj(u- V)M 



d t Sm - D V 2 Sm 



(5) 



±Mx(^.ftM)-^ 



where A is the tensor defined by A a p = jjt 5Zj(M x 
<9iM) Q (M x diM) p [13 . The right hand side of Eq. © 
acts as a source of itinerant spin dynamics inside the 
domain wall. Note that the spin diffusion (cx V 2 <5m) 
was neglected in Ref. Q- In Eq. (J5j) , by taking 
dt ~ 1/tm and V ~ 1/W (tm and W are respec- 
tively the domain wall dynamics timescale and width) , we 
have l/r cx ,l/r s f,6j/VF » 1/t MiV /(t m W 2 ),D /W 2 . 
Therefore, the left-hand side of Eq. §5§ is dominated 
by the torque and relaxation terms (oc l/r ox , l/r s f) and 
the lowest order itinerant spin density reduces to 



Sm 



where f3 
get 

1 



^--^-(1 - Amx)M x (u • V)M 

■0 2 M S K Ms ' y ' 
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TM 



V 



W 2 ' T M W 2 



(0) 



T ox/ T sf- By injecting Eq. ^ into Eq. ([5]), we 



1 



-Sm x M 



— = bj(u ■ V)M - ^d t M (7) 
bj 



-L-M x (A ■ d t M) + A 2 x -f V 2 [M x (u • V)M 



M s - ' ' ox M s 

Here X 2 X = Z?ot ox /(1 + (3 2 ) is the transverse spin dif- 
fusion length already mentioned in metallic spin-valves 
21 1. The term dt5m (oc t cx /tm) as well as the high- 
est order terms (oc /3 2 , (3rj, A 2 x /3, A 2 X ?7 << /?) have been 
disregarded since we retain only terms comparable to /3. 
Equation ([7]) can be easily manipulated to get the first 
order corrections to the spin torque exerted on the mag- 



netic texture 0], 
where 



■5m x M = T re 



M, 



■d t M 



P^M x dtM, 



T st =bj(u-V)M-f3jj-Mx(u 
+A c 2 x ^-V 2 [Mx(u-V)M] 

T d = -JU* x (A -d t M). 

Ms 



V)M 



(8) 

(9) 
(10) 



Therefore, the torque can be decomposed into three con- 
tributions: renormalization torque (oc dtM), spin torque 



3 



(oc (u • V)M) and damping torque (oc A ■ <9(M). In the 
remaining of this work, we neglect the renormalization 
torque (n s << M s ). 

In the following we study the combined influence of the 
spin diffusion and SMF-induced damping on two typical 
magnetization patterns: a Bloch wall and a free vortex 
core. The magnetization dynamics is governed by the 
Landau-Lifshitz-Gilbert (LLG) equation 



d t M = -7M x H cff - 
H c ff = H K ^e x + 2 



- — M x d t M + T d 

A M z 



T st ,(ll) 
-H, (12) 



where Hk is the uniaxial anisotropy field, A is the ex- 
change constant, and H is the external applied field. In 
a perpendicularly magnetized nanowire, a Bloch wall is 
described by 6{x) = 2 tan" 1 e z / w ', = 4>{t). The coupled 
dynamics gives 



d T X = jW 



^sin20 +(a + ^ 



-bj 



d T cj> = 7 



1 _l M , 2l 1 \ 1 A ox/ a 

rr Hk . „ ,1 . bj 

H z ~a r y sin 20 H 

'2 ^ W 



5W 2 



1 A, 



3W 2 



(13) 
14) 



2^2)), and X is the position of 



where t — t (l + a(a 
the domain wall center. At the lowest order, we obtain 
the velocity below and above the Walker breakdown 



,<WB 



,>WB 



W 1 H Z -bj(/3 



2W 2 I 



a(l + a 2 + rf) 
jW{a + i)H z - bj 



1 



if 



(15) 
(16) 



where ?/ = 2a?^/3A 2 is the reduced SMF parameter. The 
spin diffusion only significantly affects the transverse wall 
velocity below the Walker breakdown by enhancing the 
effective non-adiabaticity. 

Alternatively, in a two-dimensional magnetic stripe, an 
isolated vortex core is described by 9(x, y) = 2 tan -1 r/ro 
for r 



+ y 2 < r , 



7r/2 for ro < r < R, and 



cf> = tan - y/x + n/2, where ro (R) is the inner (outer) 
radius of the vortex core. We use Thiele's description of 
rigid domain wall motion [22|, where dtM = — (v ■ V)M. 
Thiele has shown that Eq. (jlip can be expressed in the 
form of the sum of forces exerted on the wall [22j]. By 
multiplying Eq. (|TTj) on the left by M- x and projecting 
the obtained equation on ~diM./^M 8: one obtains 



H 



off • 



d-M + G x (v + bju) + V ■ (av + (3bju) 



bj ^ x [MxV 2 (Mx(u.V)M)]^,M 



7 M S 3 
V 



[A-(v- V)M] • dM 



dn = o 



(17) 
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FIG. 1. (Color online) Normalized velocity of the Bloch wall, 
below v^ WB (a) and above Walker breakdown v^ WB (b). 
The parameters are a = 0.005, /? = 0.01 and r] = 0.2 nm 2 . 
The dots represent the normalized velocities in the presence 
of SMF-damping. 



M 

G = -sin0(V0 x V(j>) 

7 

Vij = -Ml^cfyQ^sm 2 + d l 9d J t 
7 



(18) 
(19) 



The integral J n d£l runs over the volume f2 of the mag- 
netic stripe. After some algebra, the longitudinal and 
transverse velocities of the vortex core read 



4 + Q c ff ffoff ■ 

4 + a c 2 ff 



where we define the effective damping as a e s 

the effective non-adiabaticity as f3 e g = C/3 -+ 
C = 2 



= Ca- 

28Ag x 



(20) 

1477 

and 



ln^. 

ro 

Fig. 1 displays the Bloch wall velocities as a function 
of the domain wall width W. The velocities are normal- 
ized to the case without spin diffusion or SMF-damping 
(Aox = V = 0). Below Walker breakdown [Fig. 1(a)], the 
velocity is moderately affected by the domain wall width. 
For example, when W = 5 nm, the normalized velocity 
increases by a factor of 2 for A cx = 0.8 nm (equivalent to 
NiFc with J ox w 0.5 cV, see also Ref. (H). The SMF- 
damping is almost absent in this regime (see dots in Fig. 
1). Above Walker breakdown [Fig. 1(b)], the velocity is 
simply not affected by the domain wall width (less than 
0.01 %). 

The case of vortex cores is very different. Fig. 2(a) 
and (b) display the longitudinal and transverse velocities 
of an isolated vortex core as a function of the core radius 
ro- While the longitudinal velocity is not significantly af- 
fected, the transverse velocity is dramatically enhanced 
by the presence of diffusion. For ro = 5 nm, the trans- 
verse velocity can be increased by a factor of 10 for a 
transverse spin diffusion length A cx = 0.8 nm. On the 
other hand, accounting for a moderate SMF-damping, 
rj = 0.2 nm 2 , the normalized velocity can be slightly re- 
duced. This observation has very important implications 
in the evaluation of the non-adiabatic torque in vortex 
walls. It indicates that the traditional way to extract ft 
must be reconsidered and that a more complete analyt- 
ical treatment of the velocities [i.e. Eq. (f2T))) ] needs to 
be performed. One way to extract the non-adiabaticity 
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FIG. 2. (Color online) Normalized longitudinal v x (a) and 
transverse v y (b) velocity of the isolated vortex core. The 
parameters are R = 150 nm, a = 0.005, j3 = 0.01 and r\ — 
0.2 nm 2 . The dots represent the normalized velocities in the 
presence of SMF-damping. 
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parameter is to estimate the polar angle tan -1 v y /v x ac- 
quired by the isolated vortex core after current injection 
[141] . Fig. 3(a) shows that this angle can be also dramat- 
ically enhanced at small core sizes due to the transverse 
spin diffusion. 

Up until now, the mechanism that was expected to 
mostly contribute to the non-adiabatic torque for sharp 
domain walls (W < 10 nm) has been the ballistic spin 
mistracking proposed in Refs. @, H^. This effect is non- 
local and increases with decreasing the exchange and 
domain wall width. To assess the importance of the 
spin diffusion mechanism compared to the ballistic spin 
mistracking mechanism, we numerically calculated the 
non-adiabaticity caused by the ballistic spin mistrack- 
ing based on Ref. Q. The non-adiabaticity parame- 
ter evaluated numerically is defined as the ratio between 
the non-adiabatic torque and the adiabatic torque at the 
center of the wall. As displayed in Fig. 3(b), the esti- 
mated contribution of the ballistic spin mistracking to the 
non-adiabaticity remains very limited. Since most of the 
ferromagnetic materials used in experiments are strong 
ferromagnets (A ox < 0.5 nm), this torque only has a siz- 
able influence for extremely sharp domain walls (W < 1 
nm). Therefore, for moderately sharp domain walls such 
as vortex cores, the transverse spin diffusion would give 
the most important contribution to the non-adiabaticity. 

Finally, let us discuss about the relevance of our re- 
sults to previous experimental data. Our calculations 
[see Fig. 2(b) and 3(a)] are consistent with the large 

in, 



141 for vor- 



non-adiabaticitics measured in Refs. 
tcx cores in NiFe. In these structures the transverse spin 
diffusion length is about A cx ps 0.8 nm [H, H| and the 
radius ro ~ 5 nm, which yields effective velocities about 
10 times larger than in the absence of spin diffusion, as 
mentioned above. On the other hand, the very narrow 
domain walls investigated by Burrowes et al. [l5[ are 
obtained for FePt, which is a strong ferromagnet with a 
very short transverse spin diffusion length, A cx < 0.5 nm. 
In this shown in Fig. 1(a) and Fig. 3(b), both 

spin diffusion and ballistic mistracking contributions are 
quenched. For strong ferromagnets, the itinerant electron 
spin is very rapidly aligned on the local magnetization 



FIG. 3. (Color online) (a) Polar angle of the vortex core as 
function of the core radius. The dots represent the normalized 
angle in the presence of SMF-damping. Same parameters as 
in Fig. 2; (b) Non-adiabaticity as a function of the transverse 
spin diffusion length A ex for decreasing domain wall widths, 
calculated based on Ref. Kjl. 



and all effects related to spin misalignment vanish. As 
a consequence, in such materials, non-adiabaticity might 
be dominated by contributions from the longitudinal spin 
relaxation 0- 
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